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We present a class of 2D systems which shows a counterintuitive property that contradicts a 
semi classical intuition: A 2D quantum particle "prefers" tunneling through a barrier rather than 
traveling above it. Viewing the one particle 2D system as the system of two ID particles, it is 
demonstrated that this effect occurs due to a specific symmetry of the barrier that forces excitations 
of the interparticle degree of freedom that, in turn, leads to the appearance of an effective potential 
barrier even though there is no "real" barrier. This phenomenon cannot exist in ID. 

PACS numbers: 03.65.Xp, 32.80.Qk 



I. INTRODUCTION 

Quantum tunneling has been one of the most im- 
portant problems in quantum mechanics since its foun- 
dation. The simplest problems of tunneling are one- 
dimensional, which is where our intuition on tunneling 
comes from. The extension of ID tunneling to many 
dimensions is not straightforward. There are many pe- 
culiarities that appear in many dimensional cases that do 
not exist in ID (for systematic studies of such differences 
see, e.g., Refs. [IHS])- Quite often many dimensional tun- 
neling is equated to the tunneling of complex (i.e., many 
particle) systems. 

Key aspects of the quantum mechanical tunneling of 
complex systems were analyzed by Zakhariev et al [H [5] 
in the mid-1960s; nevertheless, this problem has become 
an area of active research only in the past few decades 
(see, e.g., Refs. [6HTQ] and references therein). Tunneling 
of a diatomic molecule has been studied in Refs. |lH - 
[16]. Mechanisms of single and double proton transfer 
have been modelled by multidimensional tunneling [17]- 
[19]. Time-dependent numerical study of tunneling dy- 
namics of a two-particle quantum system with an inter- 
nal degree of freedom has been analyzed in Ref. [20 , and 
an enhancement of the tunneling probability due to the 
formation of a long-lived resonant state of the system in 
the barrier region has been discovered (similar analytical 
studies have been done in Ref. [21 ). It has also been sug- 
gested that collective tunneling of electrons may have an 
important contribution to multiple ionization of atoms in 
a superstrong laser field [22^ . Quantum tunneling of com- 
plex systems is not only of theoretical interest. Recent 
experiments where this phenomenon is observed directly 
include tunneling of a singe hydrogen atom [23 , resonant 
tunneling of Cooper pairs [24], and a bosonic Josephson 
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junction consisting of two weakly coupled Bose-Einstein 
condensates in a macroscopic double- well potential (25) . 

We present a class of 2D systems which has a coun- 
terintuitive property that contradicts the semi classical 
intuition: A 2D quantum particle "prefers" tunneling 
through to flying above a barrier. According to our anal- 
ysis, such "paradoxical" dynamics is caused by a peculiar 
symmetry of the barrier that leads to excitations of an 
interparticle degree of freedom. There is no ID counter- 
part of such systems. 

The rest of the article is organized as follows: In Sec. 
im we present the systems and describe the counterintu- 
itive effect. The observed "paradox" is explained in Sec. 
|III| Connections between the phenomenon and classical 



physics are discussed in Sec. IV Concluding remarks and 
a possible application of the effect to quantum control are 
presented in the last section. 



II. FORMULATION OF THE "PARADOX" 

Let us consider a particle moving in 2D (coordinates 
xi and X2) toward a barrier located at the origin xi = 
X2 = 0. The initial velocity of the particle is chosen to 
be directed along the diagonal xi = X2, incident on the 
barrier from the third (where xi < and X2 < 0) to the 
first quadrant (where xi > and X2 > 0); see Fig. [l] 
While the numerical calculations are done for a specific 
Hamiltonian, the analytical analysis that follows relies 
exclusively on the symmetry properties of the 2D poten- 
tial, making our conclusions, drawn from the numerical 
analysis, general. 

The model Hamiltonian for our system is chosen as 
(atomic units are used throughout) 



^nN{(^;xi,x2), (1) 



1 fd^ 

2 \dxl 

^n{(^; Xi, X2) = aV{xi) + 3V{x2) + Un{x2 - xi){2) 



8x2 



where N = 1, 2, 4 and a being an arbitrary real parame- 
ter. The potentials V{xi) and V{x2) describe the poten- 



2 



tial barriers near the origin, for the motion along each of 
the two coordinates. The parameter a ahows us to vary 
the relative height of the barriers. We have chosen 



V{x) = xexp(— x^), 



(3) 



which corresponds to a potential barrier preceded by a 
potential well. 

The potential Un{p) describes the coupling between 
the two degrees of freedom. In the absence of Un{p) 
[i.e., for Un{p) — 0], the 2D dynamics breaks into two 
uncoupled ID motions. Nontrivial features in tunneling 
appear as the result of nonzero coupling of the two de- 
grees of freedom. 

Before we describe the choice of Un{p) in our model, 
let us introduce the center of mass {R) and relative {p) 
coordinates 

R={xi^X2)/2, p = X2-xi. (4) 
The Hamiltonian ([T]) in these new coordinates reads 
-1 -1 



- ftN{(^;p, R), 



2MdR^ 2pdp^ 
nN{a]R,p) = aV{R- p/2)^3V{R^ p/2) 

^Un{p)^ (5) 

where /i = 1/2 and M = 2. 

Now we can specify the potential that couples the two 
degrees of freedom and see its role in the problem. If 
Un{p) is attractive, as it is in our calculations, it may 
support bound states. These bound states, and their 
symmetries, play a key role. 

Here, we set Un to describe a short-range attraction. 



[/iv(p) = -Aexp(V/d) 



(6) 



Varying the parameter tat, we change the number of 
bound states supported by the attracting potential. In 
the calculations, we use A = 2 and ri = 1, r2 = 1.961, 
r4 = 3.162 corresponding to one, two, and four bound 
states supported by the Hamiltonian —1/ {2p)d'^ / dp^ + 
Un{p)- The energies of these states are —0.955 for /7i; 
-1.377 and -0.372 for U2; and, finally, -1.590, -0.856, 
-0.308, and -0.012 for U4. 

Following Ref. [20], we study the tunneling within 
the time-dependent approach solving the time-dependent 
Schrodinger equation. 



id/dt — HN{a) ^Ar(Q^; a^i, X2) = 0, 



(7) 



with the initial condition at t = that reads in the coor- 
dinates Q as 

*^(a;0,i?,p) = c</.,(p)e-(«-^)'/(2-H)+^v^2^^^. (8) 

Here C is a normalization constant and is 
the ground state of the interparticle Hamiltonian, 
— l/{2p)d'^/dp^^UN{p)' In all our studies, we set m = 1, 



^=—55, (7/^ = 3, and Ecm = 1 (all values are in atomic 
units). 

The initial wave function ([8| is localized in the third 
quadrant, and we calculate the probability of finding the 
particle in the first quadrant, i.e., the probability of tun- 
nehng at later time r. The reason for using 0c/ (p) as 
the relative coordinate part of the initial wave function 
is that we wanted to avoid spreading of the wave packet 
along p before it reached the potential barrier. 

We also present the initial expectation value of energy 

En = (^Ar(a;0,Xi,X2)|^Ar(a) |^Ar(a; 0, Xi, X2)) , (9) 
El = 0.05911, E2 = -0.3631, E4 = -0.5766 

(all values are in atomic units). Rigorously speaking. 
En depends on a; however, this dependence is very weak 
because the initial wave function ([8|, independent of a, is 
mostly localized in the region where the potential barrier, 
aV{xi) + 3^(0:2), vanishes. 

The probabilities of tunneling, disintegration (see be- 
low for the clarification of this term), and refiection are 
defined as follows: 



dxi / dx2 |^Ar(a;r, Xi,X2)| (,10) 



00 

dxi 
J —00 



/O POO 
dxi / dx2 |^Ar(a;r, Xi,X2)|^ 
-00 -'O 

dx2 |^Ar(a;r,Xi,X2)|^)-l) 

) 

dxi / dx2 |^Ar(a;r, Xi,X2)|^12) 

-00 J —00 

However, since the potential barrier, aV{xi) -\-3V{x2)^ 
has "well" and "hill" regions, we also employ the corre- 
sponding "shifted" probabilities to exclude regions were 
the potential barrier is localized 

POO POO 

p[^\a,T) = dxi (ix2 |^Ar(a;r,Xi,X2)|^ (13) 
Js Js 

/ — 3 POO 
dxi / dx2 |^Ar(a;r,Xi,X2)|^ 
-00 Js 

POO P-3 ^ 

+ / dxi I dx2 |^Ar(a;r, Xi,X2)| (,14) 

r-3 



-3 P-3 

p\!^\a,T) = / dxi I dx2 |^Ar(a;r, Xi,X2)|t5-5) 
-00 J —00 



(16) 



pi^^ is the probability that a particle is trapped in the 
neighborhood of the potential barrier (in fact, mostly in 
the well region of the potential barrier). We introduce 
these quantities to verify that our conclusions are not due 
to variations of the probability density in a neighborhood 
of the potential barrier (see Ref. [26]). 

The Hamiltonian ([T]) can be viewed as the Hamiltonian 
of two ID particles, where xi^2 are coordinates of the first 
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FIG. 1: (Color online) Plots of the potential Hat (a; R, p) [Eq. 
([5|] for different a and N . Roman numerals in plot (a) label 
quadrants. Black solid lines denote the level sets of the ex- 
pectation value of the initial energy, En [Eq. ([o])], i.e., the 
boundary between the classically allowed and classically for- 
bidden regions. N — 1 and a = 3; (b) N — 1 and a — —3; 
(c) = 2 and a = 3; (d) = 2 and a = -3; (e) iV = 4 and 
a = 3; (f) = 4 and a = -3. 



probability) and does not "feel" any barrier in Fig. [TJb) 
(the particle moves solely in the classically allowed re- 
gion) while traveling along the axis p = 0, then one 
would intuitively conclude that the probability of find- 
ing the particle in the first quadrant in Fig. [TJa) ought 
to be smaller than in Fig. [ijb). By the same token, the 
very same probabilities in Figs. IT|c) andjlje) should be 
smaller than in Figs, [ijd) andjlpj, respectively. Further 
discussions of the phenomenon from the point of view of 



classical trajectories are presented in Sec. [IV 

^ are obtained from 



The results presented in Figs. [2 
the numerical solution of the time-dependent Schrodinger 
equation ^ by means of the split-operator method with 
an absorbing boundary condition. Figures |2}|4] show the 
dependence of the probabilities of tunneling, disintegra- 
tion, and reflection as function of the parameter a that 
characterizes the asymmetry of the potential barrier. Dy- 
namics of tunneling processes occurring in Figs. [TJa)-(f) 
are visualized as animations, which are available for view- 
ing in Ref. [26 . 

Remarkably, while our qualitative conclusion reached 
regarding Figs. [TJa) andjljb) is indeed correct (see Fig. 
|5|, the conclusions regarding the comparison of Figs. [TJc) 
and [TJd) and Figs. [TJe) and [TJf) turn out to be com- 
pletely wrong. In other words, the particle does prefer 
to "go" through the barrier [Figs. Illc) and file)] rather 
than flying above the barrier [Figs.^ld) andme)]. Fur- 
thermore, even though the potentials l]i(±3;i?, p) look 
similar to 1^2,4(^3; p), the particle favors motion above 
the barrier [Fig. fllb)] rather than penetration through 
the barrier [Fig. ma)] for the former pair of the poten- 
tials. This "paradox" is resolved in the next section. 



III. EXPLANATION OF THE EFFECT 



and second particles, respectively. This interpretation is 
crucial to explain the observed effect. Utilizing such a 
point of view, quantities [Eq. (11)] and p^^^ [Eq. 

(14)] can be indeed labeled as the prolDabilities of disin- 
tegration because the particles are flying apart (i.e., the 
two particle system is disintegrating) if after sufficiently 
long time r either xi > and X2 < or xi < and 

X2 > 0. 

Before stating the results of numerical calculations, let 
us qualitatively analyze dynamics of the system within a 
semiclassical consideration. Figure [l] presents the plots 
of the potentials (|5|. According to the initial condition 
[see Eqs. ([8| and ^], the particle is located on the axis 
p = and its initial velocity is directed along this axis 
toward the first quadrant, and the amplitude of the ve- 
locity is chosen such that the total energy of the particle 
equals En] hence, the boundaries between the classically 
allowed and classically forbidden regions are drawn by 
solid black lines in Fig. [l] Now compare Fig. [TJa) with 
Fig. [ijb). Since the semi classical counterpart of our 
quantum particle "experiences" the barrier in Fig. [ija) 
(penetration though a barrier is of exponentially small 



The "paradox" posed in Sec. [n]is resolved in this sec- 
tion by analyzing a perturbation theory solution of the 
Schrodinger equation. 

In this section, we shall study a general two- ID- 
particles system given by the Hamiltonian 



H 



-1 

2M dR^ ~ 
U{p) + Vi 



p 



R 



mi 



Vn R 



_P_ 
1712' 



P ,(17) 



which is already written in the (general) center of mass 
(R) and relative {p) coordinates, 

R {mixi + 7112X2) /{mi + 7712), p = X2 - xi, 
p = mim2/M, M = mi+m2, 

where Xi^2 (^1,2)7 as previously, being the coordinates 
(masses) of the first and second particles, respectively. 
Let us introduce the following notation 



U{tf,ti) = Texp -iH{tf-ti) 



(18) 



Uo{tf,U) = Texp -i{H-Vi-V2){tf-U) ,(19) 
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FIG. 2: (Color online) Probabilities of tunneling as a function 
of the height of the barrier (a), (a) P^^\a, 150) [Eq. (10)]; 
(b)pf)(a,150) [Eq. (El)]. 




(a.u.) 



FIG. 3: (Color online) Probabilities of disintegration as a 
function of the height of the barrier (a), (a) P^\a, 150) 
[Eq. In)]; (b)pf)(a,150) [Eq. 



FIG. 4: (Color online) Probabilities of reflection as a function 
of the height of the barrier (a), (a) P^\a, 150) [Eq. ([l2|]; 
{h)pi''\a,150) [Eq. (Fs}]. 



for total and unperturbed propagators, respectively. The 
sum of the potential barriers, Vi + shall be considered 
as a perturbation. The eigenfunctions, |n), and eigenval- 
ues, of the internal motion are the solutions of the 
problem 



2/i V 



\n) = E„ \n) , 



(20) 



where the index n denotes bound and continuous states. 
Introducing \nk) = |n) (g) where \k) is an eigenfunc- 
tion of the free motion of the center of mass {R\k) = 
ex.-p{ikR)/\^27r - a plane wave and |n) is an eigenstate of 
the internal motion, the unperturbed propagator reads 



Uo{tf,U) 



^ j dke~'(^-^™)^''~''^\nk){kn\ (21) 



The total propagator is a solution of the Lippmann- 
Schwinger equation written in the "post" form 

- ij^ dtUo{tf,t)[Vi^V2]e''U{t,U),{22) 
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FIG. 5: (Color online) The 



P^''\a,T) [Eq. (hol] 



and p[^\a 



of tunneling 
as functions 

of time r. (a) Comparison of Fig. [ija) vs. Fig. [ijb); (b) 
comparison of Fig. file) vs. Fig. mdj; (c) comparison of Fig. 
[Ije) vs. Fig. ^f). 



probabilities 
Eq. 

vs 



;ies c 

Si 

s. Fii 



where Wnn'{k-k') = {k'n'\ \nk) and 1^1,2) are the 

first- and second-order corrections, respectively. Higher- 
order corrections can be derived in a similar manner, but 
what is important for our further analysis is that they 
are functions of W . 

We simplify the matrix element W by representing it 
as follows 



x[Vi{q)S{R- i^p/mi - q) 
^V2{q)S{R^fip/m2-q)], (26) 



where (l)n{p) = After trivial integration over i?, 

we obtain 



Wnn-{k - k') = ^„„, (^^[fc - k']^ I ge'('=-'=')'T/i(q) 
(4^^' " ^0 / ^^''"""'^'^2(9), (27) 

where the quantity. 



where we set e ^ 0. 

Assuming that the initial condition |^(ti)) = |n) (g) 
\ipin)j where \n) is one of the eigenstates ([20]) and {i/jin) 
is a wave packet localized before the barriers [see, e.g., 
Eq. (Isl)], we obtain 



,{p) = J dpe'^''4>*Ap)4>n{p), (28) 



is called the form factor, and it is well known in the 



|^(+oo)) ^ t/o(+oo, -00) l^(-oo)) + 1^1) + 1^2) , (23) scattering theory. Its physical interpretation is the prob- 

2 ' / 2 \ ability amplitude of transferring a momentum p from the 

1^^^ = —2711 V] / dkdk'6 ( En' H En ) center of mass to the interparticle degree of freedom by 

J J \ 2M J making the transition n ^ n' . 

n' 

X |n /c ) Wnn'ik — k){k |V^in) 7 (24) Now we consider the case of identical particles: mi = 

1*2) = 4 fdkdk'dk"\n"k"){k\i^in) r iJl^'^h Then there are two 

' ' J J / \ ir / types 01 the eigenstates 01 the mternal motion: even (+), 

^n(-p) = ^n(p), and odd (-), ^n(-p) = -K{p)' Since 

XS {En- + k''^/[2M] -En- k^/[2M]) ^nn'{-p) = ^nn'{p) Wnn'{-p) = -^nn'{p)] in the CaSC 

W > nik' - k")W >(k - k') of the same (different) parity, Eq. 

" ^„ + fcv"(2M)-^„,-"fc-V(2M) + zO ' ® ^^l^^^ 



w ^ [ '^nn- {[k - k']/2) f ge'('=-'=')«[T/i(5) + 1/2(9)], if <Pn and have the same parity . 

""'^ ^ 1 ^„n' {[k - k']/2) J ||e'('=-'=')«[l^i(9) - V2{q)], if <^„ and have different parities, ^ ' 



r 



which is the product of the form factor and the Fourier transform of either the sum of the barriers or the dif- 
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ference of the barriers, depending on the parities of the 
initial and final states. 



Two conclusions can be readily drawn from Eq. (29): 



First, considering tunneling within the time-independent 
picture, Amirkhanov and Zakhariev [5 have discovered 
the violation of the barrier penetration symmetry for 
complex particles, i.e., the penetration of composite par- 
ticles through asymmetric barriers in opposite directions 
may differ. [Note that the rates of tunneling of an el- 
ementary (structureless) particle are exactly the same 
in both the directions within the time-independent ap- 
proach.] The situations when the system approaches the 
barrier from the left and from the right differ only by 
inversion of the sign of the momentum of the center of 
mass. The only part of the wave function (23) that main- 



tains the dependence on the sign of the momentum is the 
matrix element (27). Thus, the discussed phenomenon of 



tunneling asymmetry is manifested in our consideration 
as a physical consequence of the property 

Wnn'ik' -k) = Wnn'ik - k') ^ V^{q) = V2{-q). (30) 



Equation ( |3Q[ ) is not only an alternative and perhaps 
faster way of achieving the main result of Ref. [5 but 
also the generalization of their conclusion for the case of 
nonidentical barriers {Vi ^ V2). 

Second, Eq. ( [29| ) basically provides an explana- 
tion of the observed anomalies related to the potentials 
r2^(a; Xi, X2) [Eq. ([5|] pictured in Fig. [l] if we recall that 
tunneling of a 2D particle in the potential r^n(<^; ^1, ^2) is 
equivalent to collective tunneling of two equal ID parti- 
cles through the potential barriers Vi{xi) = aV{xi) and 
^(^2) = 3V(x2). Hence, Eq. (30) determines the se- 



lection rule for transitions between states of the internal 
degree of freedom induced by (collective) tunneling. The 
key point is that the probability of collective tunneling 
strongly depends on whether an excitation of the internal 
degree of freedom is possible. If a system is initially in 
the ground state and the excitations are allowed, then by 
going to an excited state, the center of mass of the system 
lowers its kinetic energy (i.e., increasing the width of the 
barrier), consequently reducing the probability of tunnel- 
ing (see Refs. [H [TTJ [121 El US] reference therein). 
We recall that the parity of the ground state is even, the 
first excited state - odd, the second excited state - even, 
etc.; therefore, according to Eq. (I29|), if 



Vi{q) = V2{q), 



then the transition from the ground state to the first 
excited state is forbidden, but if 



Vi{q) = -V2{q), 



(32) 



then the transition is allowed. Paraphrasing, we note 



that if condition (31) takes place then the interparticle 



degree of freedom may stay the same while the center of 
mass traverses the barriers, but if condition (32) holds 



change. In essence, this is the core of the observed phe- 
nomenon in Sec. ITTl 

Indeed, condition ( [32] ) is satisfied for Figs, [ijd) and 
[TJf). Hence, the tunneling probability is less in these 
cases than in cases Figs. [TJc) and[TJe) for which equality 
( [31] ) takes place. On the whole, the same conclusion is 
valid as long as the potential Un{p) can have at lest two 
bound states. In Figs. [TJa) and [ijb), when there is 
only a single bound state supported by the intraparticle 
interaction, the simple intuitive picture holds. Why is 
this the case? After all, there are also continuum states 
of the intraparticle motion available for the excitation. 

To answer this question, let us look at collective tun- 
neling from the point of view of the (time-independent) 
multichannel formalism, which is the most common 
method employed to the problem at hand (see, e.g., 
Refs. OEl El El [nHEl [27H29]). According to the 
multi-channel approach, using the expansion ^(i?, p) = 

(j)n{p)Xn{R)j the stationary Schrodinger equation, 

J n 

H^{R^ p) = E^{R^ p), is reduced to the following system 
of ordinary differential equations for unknown functions 

Xn{R)- 



1 dPXn{R) 

2M dR^ 



{R)Xn'{R) = {En-E)Xn{R), 



where Znn'{R) being effective potentials, 



dp Kip) 



Vi R 



mi ' 



R 



m2 



(t>n'{p)- 



(33) 



Such an effective potential may be interpreted as the po- 
tential barrier that the center of mass encounters while 
incident in the state n' and reflected or transmitted in 
the state n. 

Let g denote the ground as well as the single bound 
state of the potential Ui{p) [Eq. (|6|], and c a low-lying 
odd (unbound) state of the continuum spectrum, which 
is normalized to the delta function. Then, Zcg{R) = 
Zgc{K) = since the g ^ c transition is forbidden in 



Fig. fUa) ; respectively, Zgg{R) = in Fig. [Tib). The 



then the state of the interparticle degree of freedom can 



first nonzero effective potentials in Figs. [TJaj and [TJb) 
(^32) are Zgg(R) and Zcg{R), respectively, and they are plotted 
in Fig. p] Taking into account that Eg — Ei ^ Ecm = 1 
(a.u.) [Eg being the ground state energy) is the kinetic 
energy of the center of mass, we may qualitatively con- 
clude from Fig. [6] that the center of mass needs to tunnel 
through the barrier [Zgg{R)] in Fig. TJa) and flies above 
the barrier [Zcg{R)] in Fig. [lib); thus, the probability 
of finding the particle in the first quadrant in Fig. ^Co) 
prevails over the probability of tunneling in Fig. [l[a). 
Finally, we note that the first nonzero effective poten- 
tials in Figs. [TJc) and [ijd) are of the same order; the 
same statement is valid in cases Figs. [TJe) andjljf). 
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FIG. 6: (Color online) Plots of effective potentials Znn'{R) 
[Eq. (|33|]. The solid line is Zgg{R) for case (a) of Fig. [l] 



The dashed line represents Zcg (R) for case (b) of Fig. [T] 



Traces of the forced excitations, which occur in Figs. 
[TJd) andjljf), can be directly observed in obtained nu- 
merical data; these are steplike structures in P^^^(— 3,r), 

pPi-3,T), P^'\-3,t), andp('V3,r) [see Figs, gb) 
and [sjc)] and a snakelike shape of the wave function 
^2,4("^3; i?, p) that emerges from the barrier (see Ref. 
[26]). Nevertheless, the data regarding Fig. [ijb) [see 
Fig. [sja)] seems not to reveal similar jumps at first 
sight. It is due to the fact that the coupling between 
bound states [the form factor (28), more precisely] is 
bigger than the coupling of a bound state to a state of 
the continuum spectrum. Regardless of smallness, these 
transitions show up in the observation that the probabil- 
ity of disintegration in Fig. [TJa) is less than in Fig. [ijb) 

[see Fig. ||that P^^^(3,150) < P^^^(-3,150) as well as 

p^'^(3,150)<p^'^(-3,150)]. 



Concluding this section, we list pivotal factors in ex- 
plaining the "paradox" reported in Sec. |ll[ (i) This effect 
cannot exist in one dimension, it requires at least two di- 
mensions, (ii) Our explanation of the effect relies on a 
natural isomorphism between systems of one 2D particle 
and two ID particles of the same mass. The essence of the 
effect lies in possibility of tuning the potential barriers 
such that the intraparticle degree of freedom is excited, 
(iii) Dynamics of tunneling crucially depends on whether 
the intraparticle potential supports one or more bound 
states (an exact number is irrelevant for the qualitative 
description). 



IV. CLASSICAL PHYSICS AND THE 
"PARADOX" 

A classical counterpart of the quantum system at hand 
[Eq. ([T])] is a mechanical system with the Hamiltonian 

J^{pi,P2;xi,X2) = {pI +P2) /2 + l^Ar(a;xi,X2), (34) 

where pi^2 and Xi^2 are canonically conjugate variables. 
One may perform the canonical transformation to rewrite 
the Hamiltonian ( 34 ) in terms of the new canonical vari- 
ables Pr, Pp and i?, p, where the latter pair being the 
center of mass and relative coordinates [Eq. Q], 



Jf {Pr, Pp; R, p) = Pl/A + P2 + i?, p). (35) 

The connection between new and old canonical momenta 
reads: pi = Pr/2 - P^, p2 = Pr/2 + P^. 

The initial condition for the classical counterpart that 
corresponds to the initial condition ([8| is xi{0) = ^2(0) = 

^andpi(O) =pi{0) = [En - ^Nia; R, R)]^^'^ . Having 
calculated classical trajectories with this initial condi- 
tion, we observe that the classical particle does not reach 
the first quadrant in Figs. [TJb),[TJd), and[ljf); it is re- 
flected back to the third quadrant. The same conclusion 
can be reach qualitatively by calculating the force that 
acts on the classical counterpart. 



Pi, 2(^1, ^2) = -dQN{(^] xi, X2) / dxi^2' (36) 



Since 



Fi{x,x) = {2ax'^ - a)e ^, F2{x,x) = {6x'^ - 3)e ^, 

the classical particle experiences the force that deflects 
it from moving along the diagonal, Xi = X2 (which co- 
incides with the axis p = 0), and pushes it toward a 
knee-like barrier located in the second quadrant (see Fig. 
[1]); hence, the particle eventually bounces off the barrier 
back to the third quadrant. 

It is noteworthy to mention a peculiarity of numeri- 
cal calculations. We have found that it is advantageous 
to employ a (fourth-order) symplectic integrator [30^^32] 
for solving Hamilton's equations in this section due to 
the following reason: A sharp and localized shape of the 
kneelike barrier leads to an unstable motion of the clas- 
sical particle. If one employs nonsymplectic integrators 
(e.g., the Runge-Kutta methods), a very tiny time step 
must be chosen in order to properly account for the in- 
fluence of the kneelike barrier; this, in fact, often leads to 
instability of the numerical scheme for a long time prop- 
agation. A physical reason of such an instability lies in 
the fact that nonsymplectic integrators do not explicitly 
conserve energy while the symplectic integrators always 
do; hence, they give a proper long-time evolution of any 
chaotic Hamiltonian system. 

The observation of this behavior of the classical coun- 
terpart casts doubt on the quantum nature of the "para- 
dox." More precisely, is it possible that an ensemble of 
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classical particles, which corresponds (in some sense) to 
the initial wave function of the system at hand ([8| , would 
mimic the observed phenomenon? As shown below, the 
answer turns out to be negative. 

Furthermore, since it is well known that the applica- 
tion of the semi classical approximation, as a mediator 
between classical and quantum mechanics, to tunneling 
often is very fruitful in shading light on the physical na- 
ture of the studied process (see, e.g., Refs. [MllMj), the 
addressed question is important as the first step toward 
the usage of semiclassical methods for the interpretation 
of the "paradox." 

Shirokov has proposed the unified formalism for quan- 
tum and classical mechanics [35] -a reformulation of both 
the theories in terms of the same physical and mathe- 
matical concepts. Crudely speaking, this formalism is 
based on the well-known fact that observables of quan- 
tum mechanics can be converted from operators to func- 
tions (i.e., to a very similar form as in classical physics) 
by means of the Weyl representation [36]. In these terms, 
the probability density of quantum states is represented 
by the Wigner quasiprobability density distribution func- 
tion [33 • Hence, in order to construct a classical ensem- 
ble that corresponds to the quantum particle, we shall 
calculate the Wigner function W(Pr, P^; p) for the ini- 
tial condition (|8|, 

W{Pr, Pp; p) = {2^f / ^ - ^72, P - p72)(37) 

x^* (i? + R'/2, p + p'/2) e<^^^' ^^^p')dp'dR', 

where ^{R^p) = ^Ar(Q^; 0, i?, p). Since ^^(p) has no ze- 
ros and decays exponentially at infinity, we shall approx- 
imate it by a Gaussian 



(l)g{p) (X exp [-pV , 



(38) 



where we set Op = 1.5 (a.u.). Substituting Eqs. (38) and 
(|8| into Eq. p7| ), one readily obtains 

W{PR,Pp;R,p) oc exp{-(i?-i?)V4-/oV<^p (39) 

Therefore, our quantum system corresponds to an ensem- 
ble of classical particles with the Hamiltonian (35), and 



the initial state of the ensemble that corresponds to the 
initial condition ([8| can be generated by considering R^ 
p, Pr, and Pp as independent normal random variables 
with means R^ 0, \/2MEcm^ and with standard devia- 
tions CTp/V^, [V^CTi?] , [^/2cFp\ , respectively. 

Results of classical simulations of dynamics of an en- 
semble of 10^ particles are compared with ah initio quan- 
tum simulations in Fig. [7| Foremost, one may notice 
how well the classical simulations reproduced quantum 
behavior in the classically allowed regions in almost all 
the cases. One observes a qualitative agreement between 
classical and quantum results in the cases of tunneling 



[Figs. [7|^a), [tJ^c), and[7|^e)]. Nevertheless, in the cases of 
the over barrier motion [Figs. [7|^b), [7]^d), and[7|^f)], clas- 
sical mechanics gives more asymmetric probability distri- 
butions than quantum mechanics. This can be explained 
by means of a simple observation that there are the knee- 
like potential barriers in Figs. [TJb),[TJd), and[TJf), which 
force a majority of classical particles to go to the fourth 
quadrant. 

Having calculated the classical probability density dis- 
tributions, we may introduce the probability of tunnel- 
ing as well as the shifted probability of tunneling anal- 
ogously to the corresponding quantum quantities [Eqs. 
(10) and (13)]. The classical probabilities of tunneling 
in the above-barrier cases are an order of magnitude 
larger than the corresponding classical probabilities in 
the under-barrier cases. This conclusion contradicts the 
results of the quantum calculations (Fig. [5|. In other 
words, there is no "paradox" in classical physics. It is 
natural since the ensemble of classical particles should 
"prefer" going above than "tunneling through" the bar- 
rier. Hence, we have confirmed that the reported effect 
is genuinely quantum mechanical. 



V. CONCLUSIONS AND DISCUSSIONS 

In Sec. [H] we presented the 2D systems, whose po- 
tentials are plotted in Fig. [l] which hold the unexpected 
property that the probability of tunneling through a bar- 
rier is larger than the probability of flying above a barrier. 
As it was clarified in Sec. |III[ this phenomenon occurs 
due to a specific symmetry of the potential [Eq. ( [32] )] that 
forces excitations of an interparticle degree of freedom, 
thus lowering the probability of tunneling. This effect 
is overlooked by the intuitive conclusion which uses the 
language of trajectories within the quasi classical approx- 
imation that the tunneling is an "exponentially harder" 
process than flying above a barrier. First and foremost, 
we note that the quasi classical approximation, being an 
elegant and insightful approach in ID, is in fact very 
cumbersome and quite often impractical in 2D. Hence, 
in most situations of interest different modifications of 
the original quasi classical approximation that make ad- 
ditional assumptions on the wave function are employed 
(see, e.g., Ref. [9j and references therein). From this 
point of view, we conclude that a quasiclassical model 
capable of explaining the reported "paradox" must not 
only rely on the language of trajectories but also include 
the quantum transitions that are at the core of the effect. 

An important undiscussed issue is the dependence of 
the reported effect on the initial condition ([8|. If we sub- 
stitute (l)g{p) in Eq. ([8| by the wave function of the first 
excited state of the interparticle Hamiltonian in the cases 
oi N = 2 and N = A (note that Ecm must be appropri- 
ately decreased such that it would be possible to talk 
about tunneling), then one may expect that the "para- 
dox" should disappear, and one would observe a conven- 
tional situation: the probability of tunneling through the 
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FIG. 7: (Color online) Classical vs. quantum mechanics for potentials plotted in Fig. [T] Classical and quantum normalized 
probability density distributions (in the logarithmic scale) at the time moment t = 150 (a.u.). Black corresponds to the 
minimum values, whereas white corresponds to the maximum. 



barrier [Figs. [TJc) and [TJe)] would be smaller than the 
probability of flying above the barrier [Figs. [TJd) and 
[ijf)]. Indeed, since the transition from the first excited 
state to the ground state is allowed because condition 
( [32| is satisfied in Figs, [ijd) andjljf), then after making 
such a jump, the center of mass gains the energy differ- 
ence; hence, it can more easily tunnel in Figs, [ijd) and 
[ijf) than in Figs. [TJc) andjlje) where this transition is 
forbidden. 

As far as applications of the effect to quantum control 



are considered, consider a system of two neutral atoms 
that interact through the dipole-dipole interaction and 
are trapped, e.g., by a dipole trap. The magnitudes of the 
atomic dipoles depend on internal states occupied by the 
atoms. The internal states of the atoms can be changed 
for each atom independently by means of a laser with an 
appropriately tuned frequency, assuming that the atoms 
have different spectra. Performing such excitations, we 
may be able to switch between the cases where either 
condition (31) or condition (32) is valid. Hence, we may 
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allow or forbid the two atomic system to tunnel through 
the trap. 



A generalization of Eq. (29) as well as Eqs. (30), (31), 



and (32) to the case of njn ^ 3) particles is a non- 
trivial question that should be addressed in the future. 
One might expect that such a generalization of the effect 
may reveal many new varieties of the phenomenon, which 
could be interesting from the point of view of quantum 
control of tunneling of complex systems. 
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